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Abstract: A simple model for calculating the planetary 

albedo scanning angle distribution above 100 
MeV is discussed and the results are compared 
with present observations. 

It is now within the realm of technological feasibility 

to study the structure of a planet's atmosphere and 

(indirectly) its magnetosphere by means of a high-energy, 

high-resolution gamma-ray detector such as a large spark 

chamber or gas Cerenkov counter mounted on an orbiting 
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satellite (see reviews by Kraushaar and Kniffen ) . 

The bulk of the radiation at energies greater than 100 Mev, 
at the limb of the planet, originates from the decay of 
it 0 mesons produced by cosmic-rays interacting with nuclei in 
the planet's atmosphere. Thus, by scanning the planet 
from the sub-satellite point to the limb with a high- 
resolution detector one may unfold from the angular 
distribution of the gamma-ray flux, the mean total 
baryonic scale height of the atmosphere. In a planet 
possessing a magnetosphere, lower energy cosmic-rays may 
not penetrate to the atmosphere and in such regions of 
low magnetic latitude, there will be a minimum in the 

gamma-ray flux whereas in high penetration regions such as 

/ 
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magnetic poles, maxima should occur. In this paper, a 
simple model of a planet with an isothermal atmosphere is 
considered, neglecting possible geomagnetic effects. For 
purposes of some comparison with existing data, a "non- 
magnetic" earth is considered. 

The geometry of the model is shown in figure 1. The 
satellite is located at point S at an altitude a=5xl0 7 cm 
from the subsatellite point on the surface. We consider a 
cosmic-ray primary entering the atmosphere from a point P 
and colliding with an atmospheric gas nucleus at point C. 
The resulting secondary gamma-ray is emitted at C and 
travels to S where it is detected. R g is the radius of 
the earth; taken to be 6.37xl0 8 cm and h is the height of 
the collision point C above the surface. The cosmic-ray 
path makes an angle x with the vertical at C where 0<x<X M 


X M 


bciny th.6 maximum allOWSd ci.ng'lS S3 ShOWH • 


X M (h) = sin”*[R e /(R e +h)] 


(1) 


The angle that the line CS makes with the vertical at C is 

denoted by x • Thus, we may write h as a function of \ , a, 
a u 

and the scanning angle \p so that h(x o ) with ^ and a given is 

a 


h { X a 7 ^ r a) = R e [(l+^-) 

e 


sin^ 

sinx a " 


( 2 ) 
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The parameters s and t shown in figure 1 are given by 
the relations 


s = (R e +h) sin X 
and t = (R 0 +a) sin'i|> . 

We denote by y , the cosine of the angle between the 
incoming cosmic ray at C and the outgoing gamma-ray. 

It then follows that 

y = cos x a cos X +sin X a s ^- n X cos a > 
a being the angle in the plane tangent to the vertical at c. 
If we assume isotropic production of pions in the cms 

system, the angular distribution of photons in the laboratory 
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system is given by 


G(y;E p ) 


= Y c (i-e c y) 


(4) 


where y is the cms Lorentz factor and B c c the cms velocity 
c 

of the cosmic-ray-nucleus system. Here, for simplicity, 
we assume that the baryons in the nucleus act as if they are 
unbound at these energies (>lGeV) so that 


G(y;E p ) 


r r 

i • 

= u 


E /M+l 
£ 


1 .- 


E p /M+1 1 
E /M+l ^ J 


(5) 


where M is the nucleon rest mass (0.938GeV). We consider 
here an isothermal atmosphere with density 


P =P o e 


-h/H 


( 6 ) 


where £ = 1.225x10" 3 g.cm” 3 and the scale height H-10 cm. 
The incoming cosmic-ray intensity is taken to be 


KE _2 ' 2 (cm 2 .s-sr-GeV)" 1 , K=0.437 


(7) 
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The cosmic-rays are taken to have an absorption mean- 

free-path A =120g/cm 2 and the photon absorption mean-free- 
P 

path is taken to be A.y=60g/cm 2 . 

Thus, the cosmic-ray intensity at C is given by 


I c (EpjX/^/a) = KE p ” 2,i i (x;Ka) 


( 8 ) 


where i(x;^/a) = exp 


p o s f ^ -(s cscA-R ) /H 2 


° f 


e ' esc Ad A. 


The resulting photons are attenuated by a factor e 
P~t 


n 


’°1 f’ , ' Xa 

A J 

Y \ 


d<j>csc 2 <j>e~ (t cscc|>-R e )/H 


-n where 


(9) 


We assume the production cross section for tt -mesons 

4 

in the atmosphere to be of the form 


a (E ) =3. 5x10 2 5 E 5 (GeV) 

TT p P 


( 10 ) 


Then, the gamma-ray intensity from pion decay 


/ r-i A »# _ tt\ 


vr. ^./uixifcjv^ reacning uie satellite yxvtau uy 


Y 


I ( ip, a) = 4Np o Kt 


I 


t r - ^ 


esc XdX_ 

Cl CL 


TT 

| da j d 

o J-\c 


d(cosX) x 
c °sX M {' 


(ID 


X 


r 


MAX , 

dE e _(T1+ H ) E 
P P 


: G ( \x ; E ) a 0 

p' TT° 


(E ) i ( x ; ^ f a) 


C 

where E c , the lower cutoff energy is, in general, a function 
of geomagnetic latitude, but is here taken to be 2 GeV, E^^ 
for purposes of numerical calculation is taken to be lOGeV 
and N is Avagadro's constant. 
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The results obtained from a numerical evaluation of 
equation (11) are shown in the polar diagram in figure 2 
with the earth superimposed to scale. They show a strong 
peaking in intensity toward the horizon due mainly to the 
forwardly-directed gamma-rays (see equation (4)). A more 
realistic pion-production model than the cms-isotropic model 
used here would result in an even stronger "limb-brightening" 
effect which increases with energy. Above the horizon angle, 
the flux drops rapidly to zero as the atmosphere disappears; 
below the horizon angle, the atmosphere becomes optically 
thick and the flux distribution approximates a sec^-law as 
indicated by the dashed line in figure 2. 

In figure 3, the calculated results are compared with 
some of those obtained on 0S0-3 (Ref. 5) . The calculated 
distribution (dashed line) is folded in with the angular- 
response function of the 0S0-3 detector, generating the solid 

5 

line. Some of the data given by Borken are shown for comparison. 
The agreement in intensity is probably fortuitous; a higher 
average cutoff energy would have reduced the calculated flux 
and the absolute sensitivity of the detector has since been 
revised (Kraushaar et al. Ref. 6). Later data on 1(40 (Kraushaar, 
et al. Ref. 6) , show a somewhat more pronounced limb-brightening 
effect as would be expected with a higher average geomagnetic 
cutoff and a more anisotropic cms-pion distribution. However 
there is generally good agreement between the 0S0-3 
data and the calculated 1(40 given here. Also in support of 
the model given here, rough data obtained by OSO-3 on the hardness 
of the spectrum of gamma-rays at the limb indicates a it 0 -decay 
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origin. The softer spectrum obtained for \p - 0 may indicate 
an admixture of cascade radiation but it should also be 
kept in mind that backward pions produce gamma-rays with 
lower average energy than forward pions (the direction 
referring to the cms system) . Also, as Borken points out 5 , 
only multiply degraded cascade photons should contribute to 
the nadir flux and these are, in general, below the threshold 
energy of the detector. In addition, Borken points out 
that a cascade flux should have a secondary maximum at the 
nadir, contrary to the observations. The problem of the nadir 
flux does not have an easy solution and must be treated by 
more sophisticated methods. 

Unfortunately, because of the limitations in the resolu- 
tion of the 0S0-3 detector, the model presented here may not 
provide a unique fit to the observations. It can clearly be 
seen in figure 3 that the shape of the observed curve is 
resolution dominated and that the broad peak observed just 
below the horizon is caused by a smearing-out of. a narrow 
bright horizon flux (dashed line) . There is clearly room for 
improvement in both theory and observation but reasonably 
good results for the simple model presented here are encouraging. 
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FIGURE CAPTIONS: 

Figure 1. Geometry of the atmospheric gamma-ray model 
as discussed in the text. 

Figure 2. Polar diagram of calculated I(^) centered on 

satellite detector with earth below shown to scale. 

Figure 3. Comparison of calculated results with 0S0-3 
observations . 
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